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Abstract. The automated detection of solar features is a technique which is relatively underused 
but if we are to keep up with the flow of data from spacecraft such as the recently launched 
Solar Dynamics Observatory, then such techniques will be very valuable to the solar community. 
■ Automated detection techniques allow us to examine a large set of data in a consistent way and 
{J*} ' hi relatively short periods of time allowing for improved statistics to be carried out on any results 
obtained. This is particularly useful in the field of sunspot study as catalogues can be built with 
sunspots detected and tracked without any human intervention and this provides us with a 
detailed account of how various sunspot properties evolve over time. This article details the use 
r—i of the Sunspot Tracking And Recognition Algorithm (STARA) to create a sunspot catalogue. 
fV^ This catalogue is then used to analyse the magnetic fields in sunspot umbrae from 1996-2010, 
. taking in the whole of solar cycle 23. 
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1. Introduction 

To examine the magnetic fields measured in sunspots it is useful to have large datasets 
as there are vast differences between a simple sunspot surrounded by quiet sun and a 
J> ' sunspot which is in the centre of a complex active region. The large dataset was assemble d 
by using an automated sunspot detection algorithm developed bv lWatson et al. (2009). 
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image processing 



The Sunspot Tracking And Recognition Algorithm (STARA) is a quick and reliable way 
to process a large number of solar images and has been tested on images from a variety 
of sources including ground based observatories (such as Kanzelhohe Observatory, see 
http://www.kso.ac.at/sonnenbeobachtung/spot_rec_en.php for details on how the 
algorithm is being used), the MDI instrument on SOHO and the HMI instrument on the 
SDO satellite. 

The data used in this article were level 1.8 data recorded by the MDI instrument 
( Scherrer et al. 1995() and are taken from the launch of the instrument in 1996 through 



to early 2010 which covers the whole of solar cycle 23. We use both white light continuum 



observations and magnetograms which allows us to detect the sunspots and determine 
their magnetic properties at the same time. 



2. Creating the Catalogue 

To ensure that the magnetic fields of sunspots could be measured at the same time 
as their detection, the times of measurement had to be as close as possible. Due to the 
cadence of the MDI measurements, this was only the case once per day at 00:00UT giving 
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a dataset of around 5000 continuum and magnetogram images to process. Processing 
these images takes 30-40 hours on a single processor depending on the number of sunspots 
present. To extract the sunspots f rom the dat a, tec hnique s from the field of mathematical 
morphology were used (see Matheron (1975) and Serra I f|l982h for more detail). 




Figure 1. An example solar continuum image from 
MDI. The black horizontal line passes through 2 
sunspots. 



Figure 2. The 

various image 
processing steps 
involved in sunspot 
detection with the 
STARA algorithm. 
Top left: inverted 
profile from the 
horizontal line 
in Figure [1] (the 
two peaks show 
the two sunspots 
on that line), top 
right: after erosion, 
bottom left: after 
dilation, bottom 
right: after sub- 
traction of original 
profile. 
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Figure [T] and Figure [5] shows the various steps involved in detecting sunspots using the 
STARA code and the process works as follows (note that this example is given in 2D 
for simplicity looking at the two sunspots along the dark line in Figure [T] but the same 
applies to the full 3D Sun case) : 

• Invert the image so that the sunspots appear as bright peaks on a darker background. 
This is shown in the top left panel of Figure [5J 

• A 'top-hat' transform is then applied which consists of an erosion and a dilation 
(explanations of these terms can be found in Serra, 1982). 

• Subtract the profile after this transform from the original to obtain the bottom right 
panel of Figure [H 

• Apply a threshold to give locations of the sunspots. 

The full 3D case works in the same way but rather than operating on a 2D 'U' shaped 
profile, a 3D bowl shaped profile is used. Further processing is then applied to separate 
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the umbra and penumbra of the sunspots. To complete this process on a single image 
takes around 4 seconds. The sunspot locations are then superimposed on a magnetogram 
taken at almost the same time so that the magnetic fields can be recorded. This is 
repeated until all of the images have been processed. 



3. Magnetic fields in sunspots from 1996-2010 

As the magnetic fields present in sunspots had already been measured, we could easily 
look at the trends present over the length of the catalogue. 

Figure [3] shows the maximum magnetic field detected in the umbra of all sunspots 
present on a given day. It is assumed that the field in the sunspot umbrae are in the local 
vertical direction and so a cosine correction is applied to the MDI line of sight magnetic 
field. To minimise the effects of this, only sunspots with a value of /i > 0.95 were used 
(where /i is the cosine of the angle between the local solar vertical and the observers line 
of sight). Also, if there is a day with no sunspots, that day is omitted from the plot. 

We can immediately see that there is a large variation in the sunspot fields even over 
short timescales of a few weeks and that the majority of maximum umbral fields fall 
b etween 1500 and 350 gauss . 

Penn fe Livingston I (|2006h looked at this long term trend using the McMath-Pierce 



telescope on Kitt Peak and measured the magnetic field in the darkest observed part 
of sunspots from around 1993 to the present. This was done by measuring the Zeeman 
splitting present in the Fe I line to infer a magnetic field strength at the location of 
the measurement. When their whole data set is taken into account they find a trend of 
decreasing sunspot magnetic field strength of about 52 Gauss per year which is obtained 
by binning the data by year and looking at the mean of each bin along with the standard 
error on the mean. In Figure [4] we show the same treatment of the data from the STARA 
algorithm. The best fitting straight line to the data is shown with a dashed line and there 
are two sets of error bars present. The thin error bars show the standard deviation of all 
the data in that bin and the thick error b ar only takes into account data which would 
meet the criteria for being a sunspot by IPenn fe Livingston" ( 2006 ). This is because 



their data excludes pores, some of which can be as large as 10" in size and we observe 
sunspots that are of that order in size. These sunspots are the primary reason for the 
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Figure 3. Maximum 
sunspot umbra fields 
from 1996 to 2010 as 
measured by the STARA 
algorithm. 
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Figure 4. The data from Fig- 
ure [3] have been binned by 
year and the mean of each 
bin plotted. The dashed line 
is the best fit straight line 
to those points. The gradi- 
ent corresponds to a trend 
of -23.6 gauss per year. Thin 
error bars correspond to the 
standard deviation of all data 
in Figure [3] whereas thick er- 
ror bars only take into ac- 
count the data which met the 
criteria f or being a sunspot 
in th e iPenn fc Livingston I 
(|2006D article. The solid line 
is the scaled international 
sunspot number over the 
same time period shown for 
reference. 



fields in Figure [3] that lie below 1500 gauss as they correspond to small sunspots that do 
not yet have magnetic field strengths comp arable to larger, fully forme d spots. 

We observe the same decreasing trend as IPenn fc Livingston I (|2006l ) but with a shal- 
lower gradient of 23.6 gauss per year. However, it could be argued that there is a slight 
cyclic variation - unfortunately we have no cycle 22 data with which to investigate this. 
The data from the STARA algorithm also has a much larger spread making the errors 
larger in this result. However, the advantage of making measurements in this way is that 
it allows for a completely automated system that always processes data in the same way. 
From this data we cannot say for certain that a long term trend exists until data from 
the new solar cycle is obtained and a scaled plot of the international sunspot number is 
included in Figure 0] to show that the mean magnetic field is increasing and decreasing 
along with the sola r activ ity. A change of 600G over the solar cycle, as suggested by 
Penn fc Livingston"! ( 2006 ) would cause a change in t he mean umbral radiu s as a rela- 
tionship h as been shown by Kopp fc Rabin (1992) and Schad fc Penn I ( 2010l) but obser- 
vation s by IPenn fc McDonald I (|2007F could not uncover this in the data. Mathew et al 



( 20071) also suggests that the size distribution of sunspots, although constant over the 
solar cycle, could introduce a bias into small sunspot samples if the size distribution of 
spots is not calculated. 
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